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ABSTRACT 

Seven giant radio pulses were recorded from the millisecond pulsar PSR B1937+21 
during w 8.1 minutes observation by the Ooty Radio Telescope (ORT) at 326.5 MHz. 
Although sparse, these observations support most of the giant pulse behavior reported 
at higher radio frequencies (430 to 2380 MHz). Within the main component of the 
integrated profile, they are emitted only in a narrow 47 /us) window of pulse phase, 
close to its peak. This has important implications for doing super-high precision timing 
of PSR B1937+21 at low radio frequencies. 
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1 INTRODUCTION 



Recently Kinkhabwala and Thorsett (gOOOJ) (henceforth KT) 
studied giant radio pulses from the fastest known millisec- 
ond pulsar to date, PSR B1937+21, simultaneously at fre- 
quencies 2380 mega hertz (MHz), 1420 MHz and 430 MHz, 
using total intensity data from the A recibo Radio Telescope 
( ART ) (see also Sallmen & Backer (1995); Cognard et al. 
(1996)). They used coherent dedispersion to obtain a time 



resolution of w 0.1 to 0.2 micro seconds (ps). Their main con- 
clusions are (1) " at each frequency, giant pulses are emitted 
only in narrow (^ 10 /is) windows of pulse phase located « 
55-70 jj,s after the main and interpulse peaks", (2) their " 
... mean arrival phase appears stable", (3) their cumulative 
intensity distribution is roughly power law with an exponent 
of -1.8. 

This article reports the result of analyzing seven giant 
pulses from PSR B1937+21, obtained at 326.5 MHz using 
ORT; only a single linear polarization was available, and in- 
coherent dedispersion was used. Although the data were ob- 
served with a sampling interval of 102.4 [is, the pulse phases 
of giant pulses were obtained with accuracies of 15 to 20 
/xs; the instrument, the method of observation, and the sev- 
eral variations of data red uction for diverse p urpose s are 
given in Vive kana nd et al. (199£), Vivekanand (2000), and 
Vivekanand (2001). While the above pulse phase accuracies 
are two orders of magnitude lower than those of KT, they 
are sufficient for the current purpose; they are much smaller 
than the widths of giant pulses that are affected by inter- 
stellar scattering and dispersion broadening at 326.5 MHz 
(discussed later). 

Table 1 lists 5 of the total 7 data files observed during 
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DATE 



UT 



PERIOD SNR POS ERR 



1995-05-16 
1995-05-16 
1995-05-22 



1995-05-22 
1995-05-22 



22:11:28 1.55770189 
22:13:45 1.55770191 
22:26:31 1.55770736 



22:28:17 1.55770738 
22:30:06 1.55770740 



9.11 
8.43 
18.44 
17.52 
10.88 
14.20 
14.85 



0.244 
0.235 
0.774 
0.231 
0.248 
0.244 
0.248 



0.010 
0.015 
0.009 
0.010 
0.014 
0.010 
0.010 



Table 1. Columns 1 and 2 contain the date and time (UT) of 
observation; column 3 contains the period (in ms) used for folding, 
derived from the TEMPO package; column 4 contains the peak 
signal-to-noise ratio (SNR) of the giant pulse, while the last two 
columns contain the pulse phases of giant pulses and their errors 
(in units of phase within the period). 



the later half of May 1995, in which giant pulses were found. 
For technical reasons, the first 6 000 of the 681 984 samples 
(each of duration 102.4 /xs) were ignored in each file. The 
data were folded at the corresponding period (table 1); 15 
times samples were synthesized in each period, whose dura- 
tion differed from file to file, but is approximately 103.85 /is. 
In each file, consists of 44 436 useful periods of data, the five 
highest peaks were found, and giant pulses were identified as 
those whose peak signal to noise ratio (SNR) was higher than 
or close to 10.0. Then their pulse phas es an d their uncer- 
tainties were derived as in Vivekanand ( 2001 ). Then the in- 
tegrated profile of each file was cross correlated with that of 
file observed at UT 22:24:37 on 1995-05-22 (no giant pulses 
were found in this file), and the data were shifted in each 
file (using FFT technique) by the corresponding amount of 
time. 

PSR B1937+21 has a very small rotation period; so it 
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TIME [ms] 




TIME [ms] 

Figure 1. Top panel: Ten periods of data centered on the 
strongest spike observed in PSR B1937+21 at 326.5 MHz using 
ORT (file observed at UT 22:26:31 on 1995-05-22). The abscissa 
is in milli seconds, while the ordinate is in units of the signal 
to noise ratio. Bottom panel: The period containing the giant 
pulse. 



requires fast sampling rates that are not available at ORT. 
It has a high dispersion measure; so coherent dedispersion 
technique is required to maintain the required sensitivity; 
that is also not available at ORT. Therefore observing this 
pulsar at ORT is a difficult exercise, although ORT is the 
most powerful telescope operating currently at « 327 MHz. 
The top panel of fig. [j] shows ten periods of PSR B1937+21 
centered around the strongest giant pulse observed at ORT. 
The ordinate was obtained in un its of the signal to noise 
ratio as described in Vivekanand (2001). The effects of dis- 
persion smearing and interstellar scattering are evident in 
the bottom panel of fig. jjj. 



2 INTENSITY DISTRIBUTION 

These data are too meager to verify the cumulative inten- 
sity distribution reported by KT of giant pulses of PSR 
B1937+21. So only a consistency check will be done here, us- 
ing their 430 MHz observations. This frequency is appropri- 
ate for comparison with current ORT observations since in- 
terstellar scattering of giant pulses affects their 430 MHz ob- 
servations, just as it does the ORT observations, but hardly 
their 1420 MHz or 2380 MHz observations. 

The effective collecting area of ART is w 49 700 square 
meters (m 2 ), obtained from the antenna gain G value of 18.0 
Kelvin per Jansky (K/Jy) quoted by Thorsett & Stinebring 
( dl990t )) for observations of PSR B1937+21. That of ORT 
is « 7400 (m 2 ); this is ob tained from the G = 1.9 (K/Jy) 



eration of ORT, and noting that the current observations 
have been done in the total power mode, which gains a fac- 
tor i/2. This reduces by factor « 6.7 the sensitivity of ORT 
with respect to that of ART. The system temperature of 
ART i s 170 K towards PSR B1937+21 (Thorsett & Stine- 
bring ( 1990[ )); that of ORT is the receiver temperature plus 
spillover of 110 K (Selvanayagam et al ( 1993)) plus the galac- 



tic background contribution of 202 K, obtained by sc aling 
117 K (see equation 8-6 of Manchester & Taylor ( 1977| )) by 
the freque ncy sc aling factor ( 400/326. 5) 2,7 w 1.73 (see Salter 
& Brown ( |l988| ) for the exponent 2.7). This reduces ORT's 
relative sensitivity by the further factor 312/170 ~ 1.84. The 
total avail able b andwidths at the two telescopes contribute 
a factor a/10/8 w 1.12, while a factor \fl arises due to the 
availability of only a single polarization at ORT. 

At 430 MHz, KT look for giant pulses by integrating 
the flux density within a 150 /is window "located on the 
tails of the main pulse (MP) and the interpulse (IP)"; in 
this work they are located by means of their peak SNR, 
measured within a single time sample of duration 102.4 fis. 
Further, interstellar scattering broadening of giant pulses in 
KT is w 30 /is at 430 MHz (their table 1), while interstellar 
dispersion smearing is negligible due to coherent dedisper- 
sion. Thus at ART the energy of giant pulses can be as- 
sumed to be effectively spread over the detection window 
width of 150/is. On the other hand, scatter broadening at 
ORT is « 30 x (430.0/326. 5) 4 « 90.0 ^s (see Manchester 
& Taylor (1977)); and dispersion smearing within each fre- 
quency channel is ~ 165 /is. Now, scatter broadening occurs 
at each frequency within the observing bandwidth. Thus 
at ORT the energy of giant pulses can be expected to be 
spread over w V90 2 + 165 2 « 188 fjs. Thus the relative sen- 
sitivity of ORT is reduced by « 188/150 « 1.25 due to 
s mearing of gi ant pulse energy over a longer time, and by 
-\/l50.0/102.4 « 1.21 due to the reduced measurement win- 
dow. 

ORT's relative sensitivity is further affected by two un- 
certain factors. The first is a reduction by factor a, due to 
refractive interstellar scintillation, phasing of ORT , etc; a 
can be ss 2 or even more at ORT (see Vivekanand (1995)). 



(Selvanayagam et al (11993)) for the correlation mode of op- 



The second is /3, due to the relative SNR thresholds cho- 
sen in both works. The threshold is « 10 in this work; i.e., 
at ORT one has observed 7 giant pulses in ~ 8.1 minutes of 
observation above a SNR threshold of w 10. In KT this num- 
ber is not explicitly stated. However, from their fig. 7 one 
can deduce that a fraction of « 10~ 4 of their « 10 6 total pe- 
riods (KT mislabel these as "pulses" in the second last para 
of page 365, which can be confused with the giant pulses) 
have integrated fluxes greater than w 2000 [Jy— /is], at 430 
MHz, within the 150/is window. The corresponding flux den- 
sity threshold will be » 2000 [Jy-fis] / 150 [/is] w 13 [Jy]. 
Now, the noise of ART within the 150/is window would be 
170 [K] / G [K/Jy] / 2 / y/10 [MHz] X 150 H w 0.12 Jy 
(the extra factor of 2 in the denominator arises due to the 
definition of the gain G of an antenna). Thus KT would 
have observed ~ 100 gia nt pu lses above the SNR threshold 
of * 110 (Kinkhabwala ( fcOOlj) ). 

Thus the relative sensitivity of ORT is enhanced by a 
factor /3 « 110/10 « 11 due to the relative SNR thresholds 
of both works. Finally, the giant pulse flux density scales 
with observing frequency as u~ 31 (KT); this would further 
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enhance the sensitivity of ORT with respect to that of ART 
by the factor w (430.0/326.5) 3 1 » 2.35. 

This leads to an effective sensitivity of ORT that is a 
factor 6.7 x 1.84 x 1.12 x 1.41 x 1.25 x 1.21 x a/11/2.35 « 
1.14 x a less than that of ART. 

Now, KT claim that the fraction of giant pulses above 
a flux density S is proportional to S~ 1S ; actually the 
abscissa in their fig. 7 is intensity 7, which is S multi- 
plied by a fixed duration of time. So the number of giant 
pulses ORT should have recorded in half an hour would be 
« 100 x (1.14 x a)" 1 ' 8 . Therefore in the w 8.1 minutes of 
actual observation, ORT is expected to record m 21/ (a) 1 ' 8 
giant pulses. 

For a « 1, this is a factor ~ 3 higher than the num- 
ber of giant pulses actually observed at ORT; a needs to be 
« 1.8 for the numbers to match, which is not unreasonable 
at ORT. Furthermore, small changes in some of the above 
numbers can cause significant change in the final result. For 
example, assuming that the scatter broadening is 40 (in- 
stead of 30) (is further reduces the number of giant pulses 
potentially observable by « 14%. As another example, if the 
giant pulse flux density is assumed to scale with observing 
frequency as v~ 2 ' 6 , instead of as v~ SA , since such a pos- 
sibility exits in fig. 10 of KT, the number of giant pulses 
potentially observable decreases by « 28%. It is therefore 
concluded that the seven giant pulses of PSR B1937+21 ob- 
served in about 8.1 minutes of the current ORT observations 
are consistent with the cumulative intensity distribution of 
KT. 



3 PHASE DISTRIBUTION 

The effective time resolution of these data is so many times 
worse than that of KT that once again a detailed comparison 
is not possible with their work. However fig. ^ shows that 
there is broad agreement. 



3.1 Mean Arrival Phase 

The mean pulse phase of the 6 giant pulses in the main 
pulse of the integrated profile of PSR B 1937+21 (fig. |) is 
0.242±0.006, in units of the pulsar period (« 1.55770 ms 
For reference, the phase of the peak of the main pulse in fig. 
is 0.247 (discussed later). The error on the mean pulse phase 
is 0.006 x 1557.7 « 9.3/is. This is obtained by taking into 
account (1) the formal estimation error on the pha se of each 
giant pulse (column 6 of table 1; see Vivekanand (2001) for 
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details), (2) the departure of the estimated phase of each 
giant pulse from the mean of all six values (column 5 of 
table 1), and (3) reducing by V6 the final answer. Assuming 
Gaussian statistics one can conclude that most giant pulses 
(in the main pulse of PSR B1937+21 at 326.5 MHz) arrive 
within a phase window of w 9.3 x 5 ss 47/xs; this is consistent 
with the narrow window of < 10 /is obtained by KT at 
higher radio frequencies. A corresponding comment can not 
be made regarding the giant pulses in the inter pulse of PSR 
B1937+21 at 326.5 MHz, since these data contain only one 
such; it occurs « 19/iS after the inter pulse peak. 

A visual inspection of fig. ^| shows that the arrival phases 
of giant pulses in the main pulse of PSR B1937+21 at 326.5 
MHz are "stable" as claimed by KT. Once again this work 
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Figure 2. Abscissa is the phase within the integrated profile, 
in units of the pulsar period. The top panel shows the integrated 
profile (in arbitrary units) of PSR B1937+21 for the « 311052 
periods of useful data in the seven files observed. The bottom 
panel shows the pulse phases of the seven strongest giant pulses 
observed in five of those files (table 1). The error bars represent 
two standard deviations. This figure should be compared with 
fig. 4 of KT. 



can not comment on the stability of the arrival phases of 
giant pulses in the inter pulse. 



3.2 Its Separation From Pulse Peak 

The mean pulse phase of the peak of the main pulse of PSR 
B1937+21 in fig. ^, obtained by weighting the abscissa with 
the ordinate in the top panel, is 0.247 with (comparatively) 
negligible error. So it is separated from the mean pulse phase 
of giant pulses by (0.242 - 0.247) x 1557.70 « -8/xs, which 
is within measurement errors. This implies that the mean 
arrival phase of giant pulses is coincident with the mean 
pulse phase of the main component of the integrated profile 
of PSR B1937+21 at 326.5 MHz. The separation of « -8/zs 
is independent of the scatter broadening at ORT; that is 
expected to bias both the mean pulse phase of a giant pulse 
(modeled as a delta function convolved with a truncated 
exponential function), as well as the mean pulse phase of the 
main pulse of the integrated profile (modeled as a Gaussian 
convolved with a truncated exponential function), by the 
same amount, viz., the scatter broadening width. 

This is inconsistent with the ~ 55 — 70/xs separation 
noted by KT. They mention that the separation has a "lower 
limit of 49 (is, at 430 MHz" ; this also appears to be the case 
from a visual inspection of the bottom panel of their fig. 2. 
Therefore the best explanation currently is that ORT ob- 
servations, obtained using only a single linear polarization, 
might bias the mean pulse phase of the peak of the inte- 
grated profile to a greater extent than they do the mean 
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pulse phases of giant pulses. This is plausible if, for exam- 
ple, the mean position angle of the linear polarization vector 
is different in the two cases; or if it changes across the inte- 
grated profile at a rate that is different for the rate of change 
across the giant pulses, since the former is expected to be 
dominated by the non-giant pulse (or the so called normal) 
emission. However this result merits verification by future 
observations. 

This analysis is obviously not possible for the interpulse 
of PSR B1937+21 in the current data. 



4 DISCUSSION 

This work supports, at 326.5 MHz, the reported behavior of 
giant pulses of PSR B1937+21 at higher radio frequencies 
(KT), although with meager data. It is encouraging that this 
was possible with data having only a single linear polariza- 
tion, and which also suffer from severe interstellar scattering. 

KT showed that timing the giant pulses of PSR 
B1937+21 gives precision comparable to that obtained by 
timing the integrated profile, on the short time scales (fs 
hours). Whether this leads to better long term timing pre- 
cision depends upon the emission characteristics of giant 
pulses - how fast does their average profile stabilize, etc. 
(KT). This problem can be studied with actual data, as is 
currently being undertaken by KT at the higher radio fre- 
quencies. Unfortunately this can not be done at 326.5 MHz 
because of the technical limitations of ORT. It would have 
been nice if one could theoretically study the problem; but 
that would require more information on the intrinsic scatter 
of the giant pulse emission phase. 
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